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Abstract 



In analogy to Gamow vectors that are obtained from first order resonance poles 
of the S-matrix, one can also define higher order Gamow vectors which are derived 
from higher order poles of the S-matrix. An S-matrix pole of r-th order at = 
I ER—iT/2 leads to r generalized eigenvectors of order A; = 0, 1, . . . , r — 1, which are also 

Jordan vectors of degree + with generalized eigenvalue {Er — iV /2). The Gamow- 
Jordan vectors are elements of a generalized complex eigenvector expansion, whose 
form suggests the definition of a state operator (density matrix) for the microphysical 
decaying state of this higher order pole. This microphysical state is a mixture of non- 
reducible components. In spite of the fact that the k-th order Gamow-Jordan vectors 
has the polynomial time-dependence which one always associates with higher order 
poles, the microphysical state obeys a purely exponential decay law. 
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1 Introduction 



The singularities of the analytically continued S-matrix that have attracted most of the 
attention in the past are the first order poles in the second sheet. They were associated 
with resonances that decay exponentially in time. In conventional Hilbert space quantum 
theory it was not clear what those resonance "states" were, since a vector description of a 
resonance state was not possible within the framework of the Hilbert space. Higher order 
poles, in particular double poles have also been mentioned, but it has long been believed 
that they somehow lead to an additional polynomial time dependence of the decay law [Q. 
However, precise derivations were not possible due to the lack of a vector space description. 

This changed when the first order poles were associated with vectors ip^ = \/2ttV\Er — iT/2) 
in a rigged Hilbert space (RHS) 0, ^, |^, 0, called Gamow vectors. They possess all the prop- 
erties that one needs to describe decaying states or resonances: These Gamow vectors are 
eigenvectors of a self-adjoint Hamiltonian |^ with complex eigenvalues zr = Er — iT/2 (en- 
ergy and width). They evolve exponentially in time, and they have a Breit-Wigner energy 
distribution. They also obey an exact Golden Rule, which becomes the standard Golden 
Rule if one replaces ip^ with its Born approximation. The existence of these vectors allows 
us to interpret resonances as autonomous physical systems (which one cannot do in stan- 
dard quantum mechanics). It also puts quasibound states (i.e. resonances) and anti-bound 
(or virtual) states |]^ on the same footing with the bound states (eigenvectors with real 
energy), which have both a vector description and an S-matrix description. Mathematically, 
Gamow vectors are a generalization of Dirac kets (describing scattering states), i.e. they 
are also eigenkets. But whereas Dirac kets are associated with a value of the continuous 
Hilbert space spectrum of the self-adjoint Hamiltonian if, the Gamow kets are not, but have 
complex eigenvalues. 

Using the entirely different theory of finite dimensional complex matrices, decaying states 
(like the K^ — K^ system) have been phenomenologically described as eigenvectors of an effec- 
tive Hamiltonian matrix with complex eigenvalues. One usually assumes that these complex 
Hamiltonians are diagonalizable 0. However, unlike hermitean matrices which have real 
eigenvalues, non-hermitean finite dimensional matrices cannot always be diagonalized, but 
can only be brought into a Jordan canonical form [IU|. Finite dimensional matrices consist- 



ing of non-diagonalizable Jordan blocks have been mentioned in connection with resonances 
numerous times in the past [[ill , |12|, |13|, |15|, and they have been used for discussions 
of problems in nuclear [|I^ and in hadron [1^ physics. Jordan blocks have also been ob- 
tained in prototypes of mixing systems [|l^, and the appearance of so-called "irreducible" 



non-diagonalizable blocks in the density matrix has been sought after for some time in 
connection with irreversible thermodynamics and the approach to equilibrium |jl6[. That 



irreducible non-diagonalizable Jordan blocks may shed light on the idea of quantum chaos 
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has been mentioned by Brandas and Dreismann |T^. Also important for the understanding 



of quantum chaos and of the statistical properties of nuclear spectra are accidental degen- 
eracies and level crossing, which in the past had been almost exclusively restricted to stable 
systems driven by hermitean Hamiltonians [|1^. Based on a finite dimensional phenomeno- 



logical expression for the S-matrix [1^, Mondragon et al. extended these discussions to 
resonance states described by a Jordan block of rank 2. 

In the present paper we shall show that the Jordan blocks emerge naturally for the matrix 
elements zb\H\iI)~) of a self-adjoint P] Hamiltonian H between Gamow vectors |%)'''^'* 
of order k = 1, 2,...,r — 1. From the generalized basis vector expansion derived here it 
follows that these r-dimensional blocks are a truncation of the infinite dimensional exact 
theory in the RHS. 

The higher order Gamow vectors \z]^)^''^ have been derived in a recent unpublished 
preprint by Antoniou and Gadella The derivation is a generalization of the method 

by which the Gamow vector (of order k = 0) was derived from the first order poles of the 
S-matrix Starting from an r-th order pole of the S-matrix element at complex energy 
z = Zji, they derived r Gamow vectors of higher order, \Eji — iV /2~Y'^\ = 0, 1, ■ ■ ■ r — 1, 
as functionals in a rigged Hilbert space. These higher order Gamow kets are also Jordan 
vectors belonging to the eigenvalue zr. 

In the present paper we generalize the RHS theory of the Gamow vectors associated 
with first order S-matrix poles, which we call Gamow vectors of order zero, to poles of order 
r. Quasistationary states in scattering experiments (i.e. states formed if the projectile is 
temporarily captured by the target) can be shown to appear not only as first order poles, 
but as poles of any order r = 1, 2, . . . (|0). In section ^ we will start from the expression 
for the unitary S-matrix of a quasistationary state of finite order r and energy En-, given in 
reference sect. XVIII. 6, and obtain from it r Gamow vectors of order = 0, 1, . . . , r — 1 
which are also Jordan vectors of degree k + 1. After a review of the case r = 1 in section |^, we 
derive in section ^ the generalized eigenvector expansion, which contains the Gamow- Jordan 
vectors as basis vectors. With these basis vectors we can give a matrix representation of H 
and of e~*^* which contains the r-dimensional Jordan blocks. In section |^, we start from 
the pole term of the r-th order S-matrix pole and conjecture the state operator for the 
hypothetical microphysical system associated with this pole. This r-th order Gamow state 
operator consists of non-diagonalizable blocks which obey a purely exponential decay law. 
This unexpected result is in contrast to the behef that higher order poles must lead to 
an additional polynomial time dependence. 

At the present time there is little empirical evidence for the existence of these higher order 
pole "states" in nature. This is in marked contrast to the fact that first order pole states 
described by ordinary Gamow vectors have been identified in abundance, e.g. through their 
Breit-Wigner profile in scattering experiments and through their exponential decay law. 
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Now that our results have obhterated the prime empirical objection of non-exponentiality 
against the existence of higher order pole states, one can continue to look for them. The 
first step in this direction is to use these higher order state operators in the exact Golden 
Rule [|] and obtain the decay probability and the decay rate, including the line widths. We 
plan to do this in a forthcoming paper. 



2 Poles of the S-matrix and Gamow-Jordan Vectors 

Since the new (hypothetical) states are to be defined by the r-th order pole of the S-matrix, 
we consider a scattering system. The S-matrix consists of the matrix elements 



= {Q-r'^\t),Q+r{t)) = {4j-{t),<p+{t)) = (2.1) 

= f dE {tp-\E-)S{E + iO){+E\(l)+) . 

J spectrum H 

Since we are interested only in the principles here, in equation (|2.1| ) (and in subsequent 
equations) we choose to ignore all other labels of the basis vectors \E^) and \E) except 
the energy label E, which can take values on a two-sheeted Riemann surface. Nothing 
principally new will be gained if we retain the additional quantum numbers b = b2, b^, ■ ■ - bj^ 
in the basis system =^ \E, h^) = \E, 62? ^3? • • • b^)^ ^ind in place of the integral over 
the energy we would just have some additional sums (or integrals in the case that some 
of the 6's are continuous) over the quantum numbers b. For instance, if one chooses the 
angular momentum basis \E^) =^ \E,l, I3, 77^), where rj are some additional (polarization 
or) channel quantum numbers (cf. [Q, sect. XX.2, XXI. 4), then ( ^.11 ) would read in detail 

= 5^ 5^ /"/" dEdE' I E', I',, r/'-) X (2.2) 



x{-E', /', r]'\E, I, 11^) {^E, /, /3,r7|0+). 

Restricting ourselves to one initial rj = rjA and one final r]' = rjB channel (e.g., riB = i]a for 
elastic scattering) we obtain 

{-E',r,l'„r]B\E,l,h,vl) = {E',l',l'„r]B\S\E,l,h,VA) (2.3) 

= SiE' - E)Si,Mvb\\S\\va) 

where 

{VB\\S\\vA) = Sr{E) (2.4) 
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is the l-th partial S-matrix element for scattering from the channel rjA into one particular 
channel tjb (e.g., the elastic channel, ris = f]A)- If we consider the l-th partial wave of the 
rjB-th channel, then the S{E) in ( |2.1| ) is given by this matrix element S{E) = S^'^{E). E.g., 
if we consider a mass point in a potential barrier, then \E^) = \E, /, l^) is the angular 
momentum basis of the mass point and, depending on the shape and height of the barrier, 
one or several resonances can exist. Many concrete examples have been studied where one 
can see how first order resonance poles = Eji. —iTi/2 move as a function of the potential 
parameters [Q. We want to consider just one pole, and in the present paper we are mainly 



interested in a higher order pole at zr. Whether physical systems exist that are described by 
higher order poles is not clear, but a few examples of second order poles have been discussed 
in the past |] [|T|]. 

With the above simplifications to one channel rjB and one partial wave /, the notation in 
is standard in scattering theory. The standard scattering theory uses the same Hilbert 
space Ti for both the set of in-states cf)^ and the set of out- "states" ■?/'". The RHS formulation 
allows us to use two RHS's for the set {</>'*'} defined by the initial conditions and the set {ip'} 
defined by the final conditions. To explain this we subdivide the scattering experiment into 
a preparation stage and a registration stage, as explained in detail in reference [^. Fig. 1 
depicts these different stages illustrating the idealized process: The in-state 0"*" (precisely 
the state which evolves from the prepared in-state 0™ outside the interaction region where 
V = H — Hq is zero) is determined by the accelerator. The so-called out-state (or -j/;""*) 
is determined by the detector; |'?/'°"*)('?/;°"*| is therefore the observable which the detector 
registers and not a state. In the conventional formulation one describes both the 0" and 
the '0'°"* by any vectors of the Hilbert space. In reality the 0" (or 0"*") and ■0°"* (or -0") are 
subject to different initial and boundary conditions and are therefore described by different 
sets of vectors belonging to different rigged Hilbert spaces. The RHS for the Dirac kets is 
denoted by 

$C7^C$^ (2.5) 

where $ is the space of the "well-behaved" vectors (Schwartz space), and the Dirac kets 
(scattering states) \E^) and l^') are elements of The in-state vectors 0'''(t) = e*^*/^0^ 
evolve from the prepared in-state 0™(t) = (^2"*")"^ 0^(^)5 t < 0, and the out-observable vectors 
tp-[t) = e'^^/^ip- evolve into the measured out-state V^°'^*(t) = {Vt')-^ ip' {t) , t > 0.[|4| 



We denote the space of {0"^} by $_ and the space of {-0 } by $+. Then, $ = $_ + $+, 
where $_ fl 7^ 0. In place of the single rigged Hilbert space ( |2.5| ), one therefore has a 
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Figure Ic. Registration of \|/°'"(0 

Figure Id. Combination of preparation and registration 

Figure 1: The preparation- registration procedure in a scattering experiment 

pair of rigged Hilbert spaces: 

0^ G $_ C 7i C for in-states of a scattering (2.6a) 

experiment which are prepared 
by a preparation apparatus, 
ij~ e $+ C C for observables or out- "states" (2.6b) 

which are registered by a detector. 



6 



The Hilbert space Ti. in ( [^.5D , ( |2.(j| a), and (2.6b) is the same, but and $_ are two distinct 
spaces of "very well-behaved" vectors. The spaces and $_ can be defined mathematically 
in terms of the spaces of their wave functions {'^E\(j)'^) and {~E\iIj~), respectively. This is the 
realization of these abstract spaces by spaces of functions, in very much the same way as the 
Hilbert space TC is realized by the space of Lebesgue square-integrable functions L^[0,oo). 
The space $_ is realized by the space of well-behaved Hardy class functions in the lower half- 
plane of the second energy sheet of the S-matrix S{E), and the space is realized by the 
space of well-behaved Hardy class functions in the upper half-plane. Thus, the mathematical 
definition of the spaces $+ and is: 



and 



7/'-G$+ iff {E\^p"''') = {-e\^Ij-) eSnm 



G$_ iff {E\(i)''') = {-^E\(P+) eSnnl 



(2.7) 



(2i 



Here S denotes the Schwartz space and S fl Ti^. is the space of Hardy class functions from 
above/below. This mathematical property of the spaces and $_ can be shown to be a 
consequence of the arrow of time inherent in every scattering experiment ||23|| . 

Being Hardy class from below means that the analytic continuation {tp'lz') of {iIj~\E^) = 
{^E\iIj~), and the analytic continuation of ('^"-E'l^"^), and therewith also {ip~\z~){~^z\(j)'^), 

are analytic functions in the lower half-plane which vanish fast enough on the lower infinite 
semicircle. (For the precise definition, see 0, ^). The values of a Hardy class function in 



the lower half-plane are already determined by its values on the positive real axis [26|. From 
and (lOI) follows that 



{^'\E-){+E\(f)+) eSnn'L , p=i 



(2.9) 



and so are all its derivatives 



{{ip-\E-){+E\(j)+)Y"^ eSnH'L; p=l, n = 0, 1, 2, 



(2.10) 



because the derivation is continuous in S. 

With the above preparations one can derive the vectors that are associated with the r-th 
order pole of the S-matrix for any value of r, in complete analogy to the derivation of the 
vectors associated with the first order poles, r = 1.0 We shall see that there are r 
generalized vectors of order = 0, 1, . . . , r — 1 associated with an r-th order pole. We call 
these vectors the higher order Gamow vectors, or Gamow- Jordan vectors (since they also 



have the properties of Jordan vectors |10]). Their first derivation from the r-th order pole 
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Figure 2: The contours in the two sheeted Riemann surface, (a) displays the contour C~ 
that results from extending the spectrum of the Hamiltonian Sp{H) = into the lower 
half-plane of the second Riemann sheet and that yields the pole term in eq. ( p.l4 ) at the 
position zr = {Eji — iV/2). (b) displays the extension of the contour into the upper half-plane 
of the second sheet with pole at which we shall not discuss here any further; it leads to 
the growing higher order Gamow vectors. 



was given in |[T9|. Here we give an alternative derivation and discuss their properties and 



applications in the generalized basis vector expansion. 

We consider the model in which the analytically continued S-matrix S{uj) has one r-th 
order pole at the position oj = zr {zr = Er — iT/2) in the lower half-plane of the second 
sheet, (and consequently there is also one r-th order pole in the upper half-plane of the 
second sheet at u; = Zr). In this paper we will not discuss the pole at Zr. It leads to 
r growing higher order Gamow vectors and the correspondence between the growing and 
decaying vectors is just the same as for the case r = 1. The model that we discuss here can 
easily be extended to any finite number of finite order poles in the second sheet below the 
positive real axis. 

The unitary S-matrix of a quasistationary state associated with an r-th order pole at 
zr = Eji—ir/2 is represented in the lower half-plane of the second sheet by (|Q sect. XVIII. 6) 



Su(u) = 2irarctan( ,,^_^_„J 2»7M _ / UJ - Er - iT /2 \ ^^^^^^ 

^ ' \uo-{ER-tTl2) ' 



-iV 



uj - {Er - iT/2) 



(2.11) 



Here, 5r{uj) = 2ir arctan( ^^^^_^^ ) is the rapidly varying resonant part of the phase shift, and 
'y{uj) is the background phase shift, which is a slowly varying function of the complex energy 
u. We have restricted ourselves to the case that the r-th order pole is the only singularity 
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of the S-matrix. Below we will mention how this can be generalized to the case of a finite 
number of finite order poles. (Note that phenomenologically only a finite number of first 
order poles have been established, but there is no theoretical reason that would prevent other 
isolated singularities on the second sheet below the real axis.) 

For our calculations we have to write (|2.11|) in the form of a Laurent series: 



Sii{uj) 



J2 Ciiuj - zrY 

l=-r 

C-r , C-r+1 



Therefore we expand the bracket in ( p.ll| ): 



+ Co + Ciiuj-zn) + ... 



(2.12) 



Su{uj) 



I J {u-ZrY 



J=0 

1=1 



-ivy 



(2.13) 



We insert this into ( p.l|) and deform the contour of integration through the cut along the 
spectrum of H into the second sheet, as shown in fig. 2a. Then one obtains 

(^-,0+) = ^ duj {^-\uj') Suiuj) {+u\(P+) + (2.14a) 

= / ciE(7/>-|E-)5„(E)(+E|0+) + (V^-,0+)p.T. (2.14b) 
Jo 

In here, Imu < on the second sheet, and 

The first integral does not depend on the pole and may be called a "background term" . The 
contour C_ can be deformed into the negative axis of the second sheet from to — cxdh. We 
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shall set this background integral aside for the moment. For the second term on the right- 
hand side of ( p.l4|) , the higher order pole term (■?/'", 0"'")p t., we obtain using the Cauchy 



integral formulas §^ ^J^^+i = ^ where = 

(^-,0+)p.T. = Mr\u^-) j^^zrj^ {M<P^) (2.16) 



n=0 
n=0 



In here, (. . . )^'2zjj nieans the ra-th derivative with respect to lu taken at the value u = zr. 
Since the kets |a;~) are (like the Dirac kets \E~)) only defined up to an arbitrary factor or, if 
their "normalization" is already fixed, up to a phase factor we absorb the background phase 
g2«7{'^) into the kets Icj") and define new vectors 

\uj'<) = |cj-)e2^^(^) . (2.17) 

(Note that the phase is not trivial since e.g. \E~)e^^^^^'^ = \E^).) For the case that the slowly 
varying background phase jIuj) is constant, the \uj'^) are up to a totally trivial constant phase 
factor identical with \uj~)] but in general (|2.17| ) is a non-trivial gauge transformation. For the 



case of a first order resonance pole, n = r — 1 = in ( |2.16|) , the phase transformation ( p.l7|) 



is also irrelevant, because for n = no derivatives are involved in ( |2.16| ). Using the phase 



transformed vectors ( p.l7 ) we can proceed in the same way as if we were using the |cc;~) with 
'y{uj) = constant. 

Taking the derivatives, we rewrite ( p.l6|) as: 

(^-,0^)p.T. = X: (-^«— i) E (l) irizl^'^ ^-'H^zn\<P^) (2.18) 

ra=0 ^ ■ ^ k=0 ^ ^ 



In here, we denote by {ip l^;'^)'-"^ the n-th derivative of the analytic function {ip \z"'), and 
with (V^"|4)(") its value at ^ = zr. Since G SnHl, it follows that 

and {ip~\z^)^'^^ are also analytic functions in the lower half-plane of the second sheet, whose 
boundary values on the positive real axis have the property e S HHl. 

Analogously, we denote by the n-th derivative of the analytic function 

Again, ^'^^'^ z\(j)~^) is analytic in the lower half-plane with its boundary value on the real axis 
being eSnnl. 

The case r = 1 (and therefore n = and = in (|2.18|) ) is the well-known case of the 
first order pole term, which led to the definition of the ordinary Gamow vectors for Breit- 
Wigner resonance states 0, ^ |, 0] . We shall review its properties in section |^. In section ^ 
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we shall then discuss the general r-th order pole term and the generalized vectors l^;^)*-'^'', 
k = 1, 2, ... ,r — 1. These vectors we call Gamow vectors of order k or Gamow- Jordan 
vectors of degree k + 1, for reasons that will become clear in section ^. 



3 Summary of the Case r = 1 

For the case r = 1 we obtain from (p.l6| ) and (|2.18| ): 



/ + 00 
dE {^-\E-){-^E\<p+) 
■OOII 



a_i 



E-{Er- tT/2) 
= -27rm_i(^-|z^)e2*^(^«)(+2^|0+) (3.1) 
= -e2^^(-«)27rr(^-|^^)(+^R|0+). 

The integral in ( p.l| ) is obtained from the integral in ( |2.16|) by deforming the contour of 
integration into the real axis of the second sheet plus the infinite semicircle and omitting 
in ( |2.16|) the integral over the infinite semicircle in the lower half-plane of the second sheet, 



because it is zero. Eq. (|3.1| ) is a special case of the Titchmarsh theorem. The value at 
z = zji of the analytic function {^~~\z~) e^^"'^^'' defines a continuous antilinear functional 
F(ip~) = {ip~\z]^) e^^"'^^'^^ = {'4'~\z]^) over the space $+ 3 and this functional establishes 
the generalized vector 1^;^) = e^*'^^^^'' G 

We can rewrite ( p.l|) by omitting the arbitrary ip' G and write it as an equation for 
the functional \z]^) E 

^ . r dE IE-) 

loon {^^r\<P' 
1 /•+°° , 1 



= — dE \E-) i — (3.2) 

' 27ry_^^^ ' ^ {+zn\<p-^) E- {En- ^r/2) ^ ' 



/ + 00 
dE \E~) 
-OOlI 



27ri y_oo„ E -zr 

over all ip' G Or we can rewrite (|3.1| ) as an equation for the operator from $_ (prepa- 
ration) to (registration) by omitting the arbitrary ip~ G and the arbitrary 0+ G 

?■ /■+°° \E^)(+E\ 

The notation for the vectors l^;^^) derives from the Cauchy theorem: {ip~\z~^) is the value 
of the function {tlj~\uj~) at the position u = zr. The definition ( p.2| ) of the Gamow vector 
is, like ( p.l|) and ( p.3|) , just another example of the Titchmarsh theorem. From the above 
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derivation, one can see why we defined the Gamow vectors: They are the vectors associated 
with the pole term of the S-matrix element. The "normalization" of the vectors 1^;^) is a 
consequence of the "normalization" of the Dirac kets \E~), and we can define Gamow vectors 
ip'-' with arbitrary normalization and phase N{zii), 

= 1;^^) N{zr) . (3.4a) 

A normalization that we shall use here is 

= \z^) (^-e2^^(^«)) V2nr . (Hb) 

The constant phase factor —e^^"'^^^\ which we introduced in ( |3.4| b) is arbitrary and of no 
significance here, and the "normalization" factor v2vrr is also a matter of convention. (Q], 
sect. XXI.4) 

The notation 1^;^) has a further meaning: It can be shown |^ that this vector is a 
generalized eigenvector of the self-adjoint Hamiltonian H with eigenvalue zr = Efi — ir/2: 

= W-l^"") = zn {^-{tp'') , VV'-e $+ . (3.5) 

where is the conjugate operator in $^ of the operator H in $. This one writes as 

H''ip^ = ZRip^ or also H'^lz^) = zr\z^) (3.6) 

or following Dirac's notation H\E~) = E\E~) 

H'ljj'^ = zr ip^ or also H\zj^) = zr \z]^) 

if the operator H is essentially self-adjoint. If one takes the complex conjugate of ( |3.5| ) one 
obtains: 

(V^^l H \r) = (V^^^l^-) {Er + t ^) (3.7) 

which one can write in analogy to ( |3.6| ) as 

(^^|/7 = z;,(V'^| or {~zr\H = z},{-zr\. (3.8) 

It has also been shown ^ that in the RHS (2.6b) the time evolution is given by a 
semigroup operator 
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(A similar semigroup time evolution operator eZ^^^^, defined however only for t < 0, also 
exists in the RHS (|2.6| a) and has similar properties.) And it has been shown that this time 
evolution operator ( p. 91 ) acts on the Gamow vectors (or on the 1^;^) G $+) in the following 
way: 

i^-le^n^n) - (e^^V-k«) = e-^^«*e-(^/^)* {r\zn) (3-10) 
or for the complex conjugate 

Omitting the arbitrary tp" G $+, this is also written in analogy to ( |3.6D as 

e-^^*^''' = e-*^-*e-(r/2)t^G (3^12a) 

or 

(V;^|e*^* = e+*^«V(^/2)i^^G| for t > . (3.12b) 

One of the most important features of the Gamow vectors is that they are basis vectors 
of a basis system expansion. To explain this we start with the Dirac basis vector expansion 
(the Nuclear Spectral Theorem of the rigged Hilbert space) which states that 

r+oo 

(j)= dE \E+){+E\<l)+) + J2 \E,n){E^\<l)) for every G $ . (3.13) 

In here, \Em) are the discrete eigenvectors of the exact Hamiltonian H = K + V, describing 
the bound states, H \ Em) = E^ \Em), and \E~^) are the generalized eigenvectors (Dirac kets) 



of H, describing scattering states [^. The integration extends over the continuous spectrum 
of if : 0< E <oo. 



Instead of the basis vector expansion ( p.l3| ) which uses Dirac kets that correspond to the 



(continuous) spectrum of H, one can use a basis system that contains Gamow vectors, and 
one obtains the so-called "complex basis vector expansion" which states: For every 0+ G 
(a similar expansion holds also for every tp' G $+), one obtains for the case of a finite 
number of first order (resonance) poles at the positions zr., i = 1, 2, . . . , N, the following 
basis system expansion: 



P — OOJI ^ 

/ dE\E^){-^E\<j>^) - 5^|^^j27rr,e2^^(^-^)(+^«j0+) (3.14) 

•^0 ^=l 

+ J2\Em){Em\<P'') for0+G$_ 
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where 



,2j7(2fl) 



^ g^^g Gamow vectors representing decaying states. 



The first integral in (|3.14|) comes from the "background term" of equation ( |2.14|) . This 
background integral is omitted in the phenomenological theories with a complex effective 
Hamiltonian |2^. The integration in ( p.l4| ) is along the negative real axis in the second 
sheet or along an equivalent contour. The third term will be absent if there is no bound 
state \Em), which we shall assume from now on. The expansion (|3.14|) follows directly from 
( p.l4| ) for r = 1 (if one assumes that the S-matrix has no other singularities in the lower 
half-plane besides the first order poles at the positions z^^, which is a realistic assumption 
if one excludes higher order poles p8|). 

The matrix representation of H in the basis of ( |3.13| ) is given by: 







Ei)\ 












En) 






E^)j 



(E^ 

^2 



V 



En 




\ 









\E,)\ 

\E2) 



En) 
±\ 



(3.15) 



< E < +00 



for aXX ^ (^^^ = ^. In ( |3.15| ), the operator is represented by a finite or an infinite 
diagonal submatrix (for a finite or an infinite number of bound states) and a continuously 
infinite diagonal submatrix, indicated by [E), where E takes the values < E < +00. If we 
consider only the case where there are no bound states (meaning we omit the submatrix of 
the Em) then the matrix representation corresponding to the basis system expansion ( ^.13 ) 
is simply given by the diagonal continuously infinite real energy matrix: 



W\E^) 



H"" \E-- 



-^) 



E 



\E-- 



^0^ 



e $, < E < +00 



(3.16) 



On the other hand, the complex basis vector expansion ( p.l4|) (again without bound 
states) leads to a matrix representation of the self-adjoint semibounded Hamiltonian H in 



the following form: 




























^Rn) 


[{H^-\E-)J 






E-)) 



ZRi 



ZR2 












^Rn 







iE)J 

V'" G C $ - oon < E <0 





'~rJ\ 








^rJ 







(3.17) 
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The same Hamiltonian H with resonances at Zji., i = 1, 2, . . . N, can thus be represented 
either as a continuous infinite matrix ( p.l6|) in the basis of ( p.l3|) , or by ( |3.17|) in the basis 



of ( |3.14[ ). The later alternative is of more practical importance if one wants to study the 



resonance properties and if one can make {ijj \E ) small. The basis vector expansion ( p.l4| ) 



is an exact representation of (f)'^ G and the matrix representation (|3.171 ) is an exact repre- 
sentation of the self-adjoint Hamiltonian. In the phenomenological descriptions by complex 
effective Hamiltonians, one uses a truncation of ( |3.14|) and ( p. IT] ), omitting the background 



integral in ( p.l4|) and the whole continuously infinite diagonal matrix (and sometimes 



even some of the zr.) in ( p.l7| ). In this approximation one represents the Hamiltonian by 
the N X N dimensional diagonal complex submatrix in the upper left corner of ( p.l7| ). For 
example, if one considers only two resonances at z^-^ = zs , -Z/ja = '^Li one then has the 
complex energy matrix: 



(3.18) 



This truncated matrix representation is only an approximation, corresponding to the ap- 
proximation of omitting the integral in ( [3.141 ). How good this approximation is depends 
upon the particular choice of the ■?/'" (or the choice of the 0^), but it can never be exact. 

4 Higher Order Poles of the S-matrix and Gamow- 
Jordan Vectors 

We shall now discuss the possibility of extending the definition of one generalized eigenvector 
l^ij)*'"'* to generalized eigenvectors of order = 0, 1, 2, . . . ,r — 1 for an S-matrix pole of 



order r. The equations ( |2.16|) and ( |2.18| ) for the pole term are rewritten (omitting 



on the right-hand side the integral over the infinite semicircle in the lower half-plane of the 
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second sheet) as 



r— 1 ^ „ 
n=0 

n=0 ■ fc=0 ^ ^ 

^ 2vrn! - 

Since G-{E) = {ilj~\E~){'^E\(f)~^) G 5 flTii, its (n + l)-st order derivatives are also elements 
of 5 n Tii , and (|4.1|) is an application of the Titchmarsh theorem in two different versions, 
for G4E) = {^Ij-\E-){+E\(I)+) and for = {{ip-\E-){+E\(f)+)Y'^\ 

The value at z = zr of the analytic functions {ip~\z'^)^'''^ (k-th derivatives of the ana- 
lytic function {ip~\z"')) defines again a continuous antilinear functional F''(ip~) = {ip~\zJi)^'''^ 
over the space 3 ip~ . The antilinearity follows from the linearity of the differentia- 
tion [{aijj^ + (3ijj2\z)Y''^ = q;*(?/']"|z)*^^^ + P*{ijj2\z)^'^''. The continuity follows because taking 
the k-th derivative D'^ is a continuous operation with respect to the topology in the space 
S DTi.^ 3 {ijj~\E"') and because is a continuous functional F. Thus F'' = o F 

is the product of two continuous maps and therefore also continuous. The continuous func- 
tionals I^^J)'-'^^ define thus the generalized vectors l-z^)*-'^^ G A; = 0, 1, . . . , r — 1. The 
r-th order pole is therefore by ([4.1|) associated with the set of r generalized vectors 



...,|4)W,...,|4)W. (4.2) 



Of the different representations of the pole term on the right-hand side of ( [4.1|) we shall use in 
this paper only the second and third line and will come back to the integral representations 
when we discuss the Golden Rule for the higher order Gamow states. 

We insert the values ( |2.15| ) of the coefficients into ( |4.1| ) and obtain 



(r.^>.T. ^ -x:(„;,)t|!-^(«,-M2.rr.i^^))_ (4.3) 

71=0 \ / / 

= - E („;i) ^2.rx: (i) (#.-14)'" '"-"(-^hi^*) 

n=0 ^ ^ ■ fc=0 ^ ^ 
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The generalized vectors ( ^4.2| ) have all different dimensions, namely (energy)" 2"*^. If one uses 
the dimensionless variable u /T as indicated in the first line of ([4.3|) , one is led to the new 
normalization of the generalized vectors 



R 



y-"-'^ = l\zl)^'^ and ^'^-<^zn\=T^^\^z/- (4.4) 



These vectors have for all values of = 0, 1, 2, . . . ,r — 1 the same dimension (energy)" 2, 
like the Dirac kets. We have in addition introduced the factor l/k\ so that these higher 
order Gamow vectors become Jordan vectors with the standard normalization. The quantity 
{ip-\z-y-^'^^ = is the value of the functional jz"^^"^ G $^ at ip- e 

However, unlike which is the n-th derivative of (tjj'lz^) G SCiTi.'^, the 

is not the n-th derivative of (■?/'" the standard Jordan vectors \z~)^^''^ are connected 
with the "derivatives" l^;")*^^) by (|4.4| ). Therefore when we want to compare our results 
with the standard results in the theory of finite dimensional complex (non-diagonalizable) 
matrices |T0[ we need to convert from the l^;")'^'^^ to the \z^)^^''\ 



With the convention ([4.4|) we obtain from 

{rA-')p.T. = -£Qi)(-^r(2vrr)X:(^-|4^^'^^""'^^^^K|0+) (4.5) 

n=0 ^ ^ fc=0 

r-l 



n=0 ^ ^ 



where we have defined the operator 

= J2 I^R^^'^ ^""'M^^iil and W;iP^ = J2 14^^'^ ^""'H^^iil • (4.6) 

fe=0 fe=0 

Here l^p'!^'' is just an abbreviation for the right-hand side of (|4.6| ), and in section |^ we will 
discuss its interpretation. Whereas Wp^^'^^ depends also upon the background phase shifts 
through ( p.l7| ), Wp%^ is just given by the S-matrix pole. 



We now return to the complete S-matrix element ( p. 14 ) and insert the pole term 
into d^g^). 



/•-ooii 

(^,0"^) = / dE{ij-\E-)Sn{E){+E\(l>+) (4.7) 

r-l . ^ n 



n=0 ^ ^ k=0 
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Omitting the arbitrary ip G and rearranging the sums in the second term, we obtain 
the complex basis vector expansion for an arbitrary 0+ G 



-oon 



dE\E+){+E\(l)+) + 



(4.8) 



A;=0 



.(fc) 




This complex basis vector expansion is the analogue of (|3.14|) if instead of S-matrix poles 
of order one (and bound states \Em)) we have one S-matrix pole of order r. To compare (|4.8| ) 
with ( |3.14|) , we write ( p.l4|) also for the case of one S-matrix pole of order one. Then using 
the same phases as in ([4.3|) and omitting all bound states and all resonances but one, 
we obtain for ( p.l4| ) 



-oon 



(4.9) 



which agrees with what we obtain from (^4.8| ) for r = 1. Comparing ( [4.8|) with (|4.9| ) or ( p.l4|) 
we see the similarities and the differences: For a first order pole there is one generalized 
vector in the complex basis vector expansion; for an r-th order pole there are r basis vectors 
in the complex basis vector expansion. Apart from the arbitrary phase-normalization factor 
— 27rr, the coefficient of the first order Gamow vector, \zjj) = \z]^) e^^"'^^'^\ has the simple 
form {'^zii\(j)'^) which resembles the component {'^E\(f)^) of the vector along the basis 
vector \E^). In contrast, the coefficients of the higher order Gamow vectors \z]^'^^''^ are 
given by the complicated expression 



r-l 



-2-r) y: 



r 

n + 1 



(4.10) 



The difference between ( |4.8|) and (|3.14]) also foretells that the role of dyadic products like 
\En) {En\ (or also of |z^)(^z/j|), which have been prominently used for pure states, will 
probably be unimportant for states associated with higher order poles. In section ^ we will 
see that for higher order Gamow states there is no meaning to being pure. 

Since the general expressions ( |4.8| ) and ( [4.7| ) are not very transparent, we want to spe- 
cialize them now to the case of a double pole, r = 2: 

-oon 



-) + 



y^'hnT{2 - t ^'^^^zn\<P^)) (4.11a) 
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If as a generalization of ( ^.4| b), we define the differently normalized Gamow vectors 
then the basis vector expansion ( [4.11| a) for the case r = 2 reads 

-ooii 



(4.12a) 



dE\E+){+E\(l)+) + 
Note that according to ([4.12| a) and ( p.l7| ) we have 



(4.11b) 



^G(i) ^ r(|2^)(i) + |z^)2^7'(2^))e2^^^(^«)v^ 



(iHb) 



and only for constant background phase shift j(^){z) = 0, n = 1, 2, ... , V^^(^) (or tp^^^^) 
given by \z]^)^^^ (or l^^lj)*^'^^). One can insert ( f4.12| b) into ( |4.11| b) and expand 0"'" in terms of 
the basis vectors and \zj^)^^^] and the same procedure one can repeat for arbitrary k to 
1) in terms of 



express 



m 



(0) 



(4.13) 



Whether the phase convention in the definition ( [4.12|) will turn out to be convenient cannot 
be said at this stage. 

The basis vector expansion can be generalized in a straightforward way to the case of 
an arbitrary finite number of poles at the positions Zf^., i = 1,2, ... ,N of arbitrary finite 
order r^. in the same way as it was done in ( p.l4| ) for = 1. This complex generalized 
basis vector expansion is the most important result of our irreversible quantum theory (as 
is the Dirac basis vector expansion for reversible quantum mechanics). It shows that the 
generalized vectors ( |4.2| ) (functionals over the space $+) are part of a basis system for the 
0+ G and form together with the kets \E~^), — oon < £" < 0, a complete basis system. 
The vectors (|4.2| ) span a linear subspace M.zr C of dimension r: 



M 



r~l 



^ = 5Zl-2K)^'^Cfc, c, G(E ^ c<i>: 



(4.14) 



fc=0 
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If there are poles at z^. of order r^, then for every pole there is a linear subspace M-zr. C 
Since the generalization to poles of order rj at energy z^i- is straightforward, we 
continue our discussions for the case of one pole of order r. 

Note that by the procedure described in this section a new label k was introduced for the 
basis vectors in the expansion ( |4.8|) , \z]i'r~-^^'' = {zn, 62, h, . . . , hJj-'r~-^^\ Usually basis vector 
labels are quantum numbers associated with eigenvalues of a complete system of commuting 
observables. That means that if in addition to H there are the M—1 operators B2, B3, . . . , Bj^ 
with eigenvalues {62, &3, • • • , bj^} = {b} = spectrum(_B2, -B3, • • • , -Bat), then the Dirac kets are 
labelled by \E,b~) = \E,b2,b^, . . . ,bjf) and in addition to the sum and integral in ( ^■13[ ) 



and ( [4 .81) , there is a sum and/or an integral over all the values of the degeneracy quantum 
numbers 62,^3, • • • , ^at, which we suppress here for the sake of simplicity. The label k of 
the higher order Gamow vectors \zR^b~)^^\ which has appeared in ( [4 .81) , is not associated 
with a conventional quantum number and is there in addition to the labels b connected with 
the eigenvalues of the set of commuting observables B2, B^, . . . , Bj^. The quantum numbers 
zr, b2, b^, . . . ,bj^ can be observed and have an experimentally defined physical meaning. It 
is not clear that the label k will have a similar physical interpretation. This means that 
(if a higher order S-matrix pole has at all a physical meaning) the different vectors l^;^^)'''^^ 
in the subspace AizR have no separate physical meaning (unless k can be given a physical 
interpretation). 

Now that (|4.^) has established the generalized vectors ([4. 2]) or the generalized vec- 



tors ( |4.13| ) as members of a basis system (together with the \E^); > E > — oon) in 



we can obtain the action of the operator H by the action of the operator on these 
basis vectors; and we can write the operator H in terms of its matrix elements with these 
basis vectors. This can also be done in the same way for any of the operators f*{H), where 
f{z) is any holomorphic function such that 

f*{H) : $_|_ — > is a -continuous operator, (4-15) 

(e.g., f*{H) = e*-'^* , f{H^) = e^^^""^ for the real parameter t > only, since for t < 
f*{H) = e*^* is not a continuous operator in $+.) For this purpose we replace the arbitrary 
G in (^73|) by ip' = f*{H)'ip~ which is again an element of $+, because f*{H) is a 
continuous operator in (by assumption ( |4.15|) ). Then we obtain by comparing powers 
of T: 

d" 

[■^j{r{H)iJ-\z^r'^--'Czn\r} = ■ 

k=0 



n = 0,l,...,r- 1 . (4.16) 
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where we have used ( |2.17| ) and 

irmriuj) = {r\f{H^)\oj-) = i^^^) (4.17) 

which follows from ( |4.15|) . The function 

G{z) = f{z){tlj-\z-){-^z\(j)-^)e'''^^'^ (4.18) 

is an element of 5 fl H^, since E S D and e^^"'^^^ as well as f{z) are 

holomorphic. Therefore we can take the derivatives G^z)^"^^ of any order 

G'(^)^"^ = E it) {fi^){r\^')f^ . (4.19) 

Inserting this into ( |4.16| ) we obtain: 

fc=0 ^ ^ 

n = 0,l,...,r- 1 (4.20) 

Since this has to hold for every 0+ G (i.e., for every G 5 fl 7Y^), it follows that 

the coefficients of each derivative '■"^'^-'(^-^rI^^) = f^^n-l) {'^ A'P^)\z=zji i^^st vanish. Thus, 

(ri/(i/^)i4>'"= (/(^)»-k^>)l, l:d*an2.-V*;':" 

By a similar argument, just comparing the coefficients of (e^*'^'-^-' (^2;|(/)+))i=^^'' rather than 
of ((^^10^))!=^^, one can show that the same equation holds for the \z]^)^^^ (with any nice 
function for 7(2;)): 



/ \ (^) 



This permits us to calculate the action of f{H^) on the generalized vectors \z]^)^''^ G 
for every f{H^) that fulfills the condition (|4.15| ). The same calculation applies to the 
generalized vectors l-^^)*-'^-' due to ( [4.21| ). Therefore we write the following equations for 
\zJi)^''^ though the same holds for \zj^)^''\ 

We first choose f{H^) = H^; then we obtain 

^ iriH^'iziy'^ = ZRirK)^'^ + Cl) irKf-'^ (4.23) 
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which can also be written as a functional equation over as 

If we use the normalization of the basis vectors defined in 
detail then we obtain 



H><\z-^y' 
H-\z-^y 



(0) 



(1) 



H-\z-^y 



Zr\Zj^?~ +i \Zji?~ 



r — 1 . 



(4.24) 



I, and write ( [4 .241 ) out in 



(4.25) 



H>^\z-^y^'^'^ = zn\z,y^'-'^ + r\z',y^-'^ 



(and the same for l^;^)^*-*^''). This means that restricted to the subspace Aizj^ is a Jordan 
operator of degree r (in the standard notation the operator is the Jordan operator of 

degree r), and the vectors \z]^)^^''\ k = 0, 1, 2, . . . , r — 1 are Jordan vectors of degree 
k + 1. |TD[ They fulfill the generalized eigenvector equation PD 



- ZR^+'lz-^y'-'^ 







We write the equations ( [4.25| ) again in the form ( [4.23|) and arrange them as a matrix 
equation. Since the basis system includes, according to ([4.8|) , in addition to the \z]^)^^^\ 
k = 0, 1, 2, . . . , r — 1, also the \E~), — oou < E < 0, we indicate this by a continuously 
infinite diagonal matrix equation which we write as: 

(^W\E-)^ = (^{i;-\H\E-)) = (i?) [{^-\E-)^ (4.26) 

where {{ip~\E~)) indicates a continuously infinite column matrix. Then ( 4.25| ) and ( 4.26 ) 
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together can be written in analogy to ( |3.17| ) as: 



{Hi; 



.(r-1) 



.(r-l) 



V {^p-\H-\E-) ) 



(4.27) 



(zr 

r 





Vo 








Zr 




\ (V-" 







(0) 

^(1) 



In this matrix representation of , the upper left r x r submatrix associated with the 
complex eigenvalue zr is a (lower) Jordan block of degree r. We have chosen the Jordan 
vectors with the normalization of ( |4.4| ) in order to obtain the Jordan block in a form closest 
to the standard form, but with F's in place of I's on the subdiagonal. 

It is instructive also to write down the adjoint (i.e. transposed complex conjugate) of 
the matrix equation (|4.27| ) because it will clarify the notation and display the upper Jordan 
block. Taking the transposed and complex conjugate of ( [4.27| ) we obtain: 



(^^''>-<~ZR\H\i;-),.../^-'^-<-ZR\H\ij-), {-E\H\r) 



/4 





r 









r 

4 








(4.28) 




\ 













{E) ) 



With the derivation of ( [4.8|) and ( [4. 271) we have reduced the problem of finding the vectors 
(and their properties) associated with the higher order poles of the S-matrix to the spectral 
theory of finite dimensional (non-normal) complex matrices, which is well documented in 
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the mathematical hterature \TU\. If in addition to the r-th order pole at zji there are other 
r j-th order poles at zj^. , then for each of these poles we have to add another Jordan block of 
degree to the matrix in ( [4.27[ ) . 

We could now refer for further results to the mathematics literature of r x r complex matrices, 
but we can also obtain these results easily from ( [4.21[) and (|4.22|) . 



Applying to the right-hand side of ( 4.21 ) the Leibniz rule we obtain 

k 



u=0 



(k-u) 



z=zr 



(4.29) 



where f^'^\z) is the z^-th derivative of the holomorphic function f(z) with respect to z and 
{ip~\z'^)'-^~'^'> = {{iIj~\z^))^''~'^^ is the (A; — i/)-th derivative of {ip~\z''). We now insert ( [4 .41 ) on 
both sides of ( [4.29|) and obtain: 



E 

u=0 



k\ 



'^k-u 



(4.30) 



From this we obtain 

k 

or as a functional equation: 



iy=0 



(4.31) 



i/=0 



(4.32) 



(Note in this calculation that {ip \z J^^*-""^ is not the n-th derivative of {ip \z )^^^\ whereas 
{ip'lz')^"^^ is the n-th derivative of {^~\z^) E S (1 1-L\. Therefore it is better to work with 
the \z-Y^^^ than with the \z->~^^\) 

In the theory of finite dimensional Jordan operators [|T^], the equality ( [4.32| ) is often called 
the Lagrange-Sylvester formula and is written as a matrix equation (using lower Jordan 
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blocks for i/^ as in (|]^): 



(4.33) 



\ 





The r X r submatrix equation of ([4.27|) is a special case of this for f{H^) = H^. Equa- 
tion ([4.33|) is not the complete matrix representation of f{H^), because the infinite diagonal 
submatrix due to the first term in 



(4.34) 



has been omitted. Equation ([4.33|) gives the restriction of f{H^) to the r-dimensional 
subspace Ai^R C. $+. 

The function of that we are particularly interested in is the time evolution operator 
can be defined in only for those values of the parameter t for which 
^iHt . — ^ -g r$^-continuous operator. This is the case for t > 0, but not for t <0. 
(For {ip^lE-) e Snn'i, the function {e'^^ip-\E-) = e-'^\ip-\E~) is an element of .SnT^^ 
only for t > 0.) Thus, for t > 0, we can use ( frj^ with f{z) = e"^^* and f^^^z) = {-itye-"\ 
and we obtain the following functional equation in M-z^ C 



\ZRr 



V I ^7 ^(*:-!^) 
R~ 



In terms of the vectors l-zj)'-^-* this can be written (using (|4.4| )): 

14) 



(4.35) 





{'^ty 1 







(4.36a) 



i/=0 
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or taking the complex conjugate (in analogy to going from ( p.l2| a) to ( ^.12| b)): 



'ZR\e 



u=0 



V (k—v) /■y 



Zr\ ■ 



(iHb) 



The vectors l^;^)'-^-' in the above equations can be replaced by the vectors \zj^)^''\ 

It is important to note that the time evolution e~*^^* transforms between different 

— V (k) 

Iz^r- , k = 1, 2... ,n, that belong to the same pole of order r at 2; = zr, but the 
time evolution does not transform out of AizR- On the basis vectors \E^) of the first term 
in (|4.8|) the time evolution is diagonal 



iEt\ 



(4.37) 



The equation ( [4.35|) and ( [4.36|) can be written as a matrix equation on the subspace M.zr C $_ 



I 6 



-iZRt 



i-itr) iz^t 

1! ^ 

01 



(-itr) izfit 

1 1 ^ 



i-itry- 

\ (r-l)! 



-izRt i-itry- 

(r-2)! 



-izjit 



(4.38) 



\ 



As an example let us consider the special case of a double pole, r = 2, /c = 0, 1. The 
formula ( [4.35| ) for the zeroth order Gamow vector is then 



-iH>'t 



and for the first order Gamow vector it is 



JO) 



(4.39) 



(4.40) 



It has been known for a long time that a double pole and in general all higher order S- 
matrix poles lead to a polynomial time dependence in addition to the exponential. However, 
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it was not clear what the vectors were that have such a time evolution. Here we have seen that 
they are Jordan vectors of degree r or less, and that they are Gamow vectors, \z]^)^^''^ G $+. 
We have also shown that the time evolution operator is not diagonal in the basis ( [4 .21 ) but 
transforms a Gamow- Jordan vector of degree {k + 1) into a superposition (|4.35|) of Gamow- 
Jordan vectors of the same and all lower degrees with a time dependence in addition to 
the exponential that depends upon the degrees of the resulting Gamow- Jordan vectors. 



5 Possible Physical Interpretations of the Gamow- Jordan 
Vectors 

In the previous section we defined the higher order Gamow vectors l^^^)'-'^-' G from the 
r-th order pole term of a unitary S-matrix. We showed that they are the discrete members of 
a complete basis system for the vectors 0^ G (|4.8|) , and we derived their mathematical 
properties: In ( [4. 24)) and (|4.25| ), we showed that they are Jordan vectors of degree k + 1, and 



in ( |4.33|) and ( |4.38|) , we obtained the Lagrange-Sylvester formula and the time evolution. 

The mathematical procedure that we used for the r-th order pole term is a straightforward 
generalization of the definitions and derivations that had been used for an ordinary, zeroth 
order Gamow vector and first order poles of the S-matrix [^]. 

Gamow states of zeroth order with their empirically well established properties (expo- 
nential time evolution, Breit-Wigner energy distribution) have been abundantly observed in 
nature as resonances and decaying states. Theoretically, there is no reason why the other 
quasistationary states (i.e. states that also cause large time delay in a scattering process 
(0 sect. XVIII. 6) and are associated with integers r > 1 in (|2.11|) ) should not exist. How- 
ever, no such quasistationary states have so far been established empirically. One argument 
against their existence was that the polynomial time dependence, that was always vaguely 
associated with higher order poles 0, has not been observed for quasistationary states. 

The question that we want to discuss in this section is, whether there is an analogous 
physical interpretation for the higher order Gamow vectors as for the ordinary Gamow 
vectors, namely as states which decay (for t > 0) or grow (for t < 0) in one prefered 
direction of time ("arrow of time") and obey the exponential law. Since we have now well 
defined vectors associated with an r-th order pole, we can attempt to define physical states 
which have well defined properties that can be tested experimentally. 

In this section we are dealing with physical questions about hypothetical objects associ- 
ated with the r-th order pole. We therefore have first to conjecture the higher order Gamow 
state, before we can derive their properties. We start with the known cases. 

In von Neumann's definition of a pure stationary state one uses a dyadic product W = 
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I/) (/I of energy eigenvectors |/) in Hilbert space. In analogy to this, microphysical Gamow 
states connected with first order poles have been defined as dyadic products of zeroth order 
Gamow vectors p9| : 

= = \zJ,){~-Zr\ = (5.1) 

(Since for the generalized vectors = \/27rr \z^) or \z]^) we cannot talk of normalization 
in the ordinary sense, it is not important at this stage whether or not to use the "normal- 
ization" factor of 27rr in VV-^. 

The time evolution of the Gamow state ( [5.1| ) is then given according to ( p.l2| ) by: 

= e-^^«* |V^^)(V^^|e^^«* (5.2) 
^ g-i(£;fl-i(r/2))t i^G^ ^^G| gi(Sfl+i(r/2))t 

= e-^'W^{0), t>0. 

Mathematically, the equation ( ^^ is to be understood as a functional equation like ( p.lO|) 
and (imi ): 

(7/>r|iy^(t)|V'2-) = e-r*(^r|iy^|^2-) or (5.3a) 
{'ip-\W^{t)\ip-) = e-^'{ij-\W^\ij-) (5.3b) 
for all ip' , ipi , ip2 £ ^+ and t>0. 

The mathematical form ( |5.3|) of the time evolution of W'^ shows how important it is in 
our RHS formulation to know what question one wants to ask about a Gamow state when 
one makes the hypothesis ( |5.1| ). The vectors tp~ G represent observables defined by the 
detector (registration apparatus). The operator W'^ represents the microsystem that affects 
the detector. Therefore the quantity {ijj'lW^lijj') is the answer to the question: What is 
the probability that the microsystem affects the detector? 

If the detector is triggered at a later time t, i.e. when the observable has been time 
translated 

— e^^*|r)(^-|e-^^* = \ij-{t)){^-{t)\ (5.4) 

then the same question for t > has the answer: The probability that the microsystem 
affects the detector at t > is 

= (^-|e-'^'*iy^e*^*|^-) (5.5) 
= e-^'{tlj-\W^\^-) . 
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This means that ( p.5| ) is the probabihty to observe the decaying microstate at time t relative 
to the probabihty {ip~\W^\ilj~) at t = 0, (which one can "normahze" to unity by choosing 
an appropriate factor on the right-hand side of ( |5.1| )). 

The question that one asks in the scattering experiment of fig. 1 is different. There the 
pole term (P.T.) of (|3.1|) describes how the microsystem propagates the effect which the 
preparation apparatus (accelerator, described by the state (p'^) causes on the registration 
apparatus (detector, described by the observable ip~). 

In conventional orthodox quantum theory one only deals with ensembles and with observ- 
ables measured on ensembles. Their mathematical representations, e.g., \4>){4>\ for the state 
of the ensemble and for the observable, are from the same space $, i.e. (p, ip & ^. 

(And if one is mathematically precise then one chooses for $ the Hilbert space, $ = 7^.) On 
this level, one cannot talk of single microsystems, and there are no mathematical objects in 
orthodox quantum mechanics to describe a single microsystem. Still, it is intuitively attrac- 
tive to imagine that the effect by which the preparation apparatus acts on the registration 
apparatus is carried by single physical entities, the microphysical systems . 



According to the physical interpretation of the RHS formulation, "real" physical enti- 
ties connected with an experimental apparatus, like the states defined by the preparation 
apparatus or the property ip defined by the registration apparatus, are assumed to be ele- 
ments of $, but states and observables are distinct. In particular, states and observables 
of a scattering experiment are distinct and described by $_ of (p.6|a) and $+ of (p.6|b). 
However, mathematical entities describing microphysical systems are not assumed to be in 
$. The energy distribution for a microphysical system does not have to be a well-behaved 
(continuous, smooth, rapidly decreasing) function of the physical values of the energy E, 
like the functions {E\ip) describing the energy resolution of the detector, or the functions 
{E\(p) describing the energy distribution of the beam. Hence, for the hypothetical entities 
connected with microphysical systems, like Dirac's "scattering states" |p) or Gamow's "de- 
caying states" \E - iT/2), the RHS formulation uses elements of $+, and [||, |9l. 
The time evolution of the "state" vectors for the decaying microphysical systems, e.g. ( ^.12[ ) 
or (|5.2|) , can be obtained from the well established time evolution of the quantum mechanical 
observable ( p.4| ) using the definition of the conjugate operator as in ( p.lO| ). 

Because of the difference between ip~ G $+ for the observables and 0+ G for the 
prepared states one needs a different mathematical description for the same microphysical 
state, depending upon the question one is asking. If one asks the question with what prob- 
ability the microphysical state affects the detector ip~{t), then the microphysical state is 
described by (|5.1|) . In a resonance scattering experiment of fig. 1 one asks another question: 
What is the probability to observe ip~{t) in a microphysical resonance state of a scattering 
experiment with the prepared in-state (p^7 

In distinction to a decay experiment, where one just asks for the probability of tp^ G 
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in the resonance scattering experiment one asks for the probabihty that relates G to 
(f)'^ G $_ via the microphysical resonance state. Therefore the mathematical quantity that 
describes the microphysical resonance state in a scattering experiment cannot be given by 
\^r){~^r\ of (|5TT|), but must be given by something like \z]^){^zii\. 

The probability to observe ip~ in the prepared state independently of how the effect of 
0"^ is carried to the detector ip' is given by the S-matrix element ( |2.14 ), , . The 
probability amplitude that this effect is carried by the microphysical resonance state is then 
given by the pole term (t/^", 0"'")p.t., equation (|2.18| ). 

In analogy to ( ^.51) one can now also compare these probabilities at different times. For 
this purpose one translates the observable ip~ in the pole term in time by an amount 
t > 0, 

ij- — ^ij'{t) = e'"'ilj-; t>0 (5.6) 

(which corresponds to turning on the detector at a time t > later than for ip' ). One 
obtains 

(V^-(t),0+)p^^^ = -27rr(e^^V-|^/j)(+^i?j0+) 

= -2vrr(V^-|e-^^*|z^)(+^H|0+) (5.7) 
= -27rre-^-*(^-|^^)(+2«|0+) 



iEiit-rt/2 



e ^ e 



P.T. 



This means that the time dependent probability, due to the first order pole term, to measure 
the observable ip~{t) in the state 0"'" is given by the exponential law: 

|(e^^V-,0+)p.T.P = e-^*|(V^-,0+)p.T.p. (5.8) 

This is as one would expect it if the action of the preparation apparatus on the registration 
apparatus is carried by an exponentially decaying microsystem (resonance) described by a 
Gamow vector. 

Therewith we have seen that there are two ways in which a resonance can appear in 
experiments and therefore there should be two forms of representing the decaying Gamow 
state (for the case r = 1 so far)[] 

by \z]^){~^zr\ in a scattering experiment, (5.9a) 
and by \z]^){~zr\ in a decay experiment. (5.9b) 



^An analogous statement holds for the Gamow states associated with the pole in the upper half-plane. 
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The first representation is the one used in the S-matrix when one calculates the cross sec- 
tion; the second representation is the one used when one calculates the Golden Rule (decay 
rate). In contrast to von Neumann's formulation where a given state (representing an en- 
semble prepared by the preparation apparatus) is always described by one and the same 
density operator W, the representation of the microphysical state in the RHS formulation 
depends upon the question one asks, i.e. upon the kind of experiment which one wants to 
perform. That a theory of the microsystems must include the methods of the experiments 
has previously been emphasized in ^0 . 



After this preparation we are now ready to conjecture the mathematical representation 
of a higher order Gamow state (a quasistationary state with r > 1). 

In analogy to the correspondence between ( ^.9| a) and (5.9b) we conjecture that for the 
case of general r we have also two distinct representations of the Gamow state. The one for 
resonance scattering is already determined as in the case for r = 1 by the (negative of the) 
pole term ( [4. 51 ), and is therefore given by 



n=0 ^ ^ ■ k=0 ^ ^ 

n— n V / 



P.T. 



n=0 

where we have used the operator defined in ( |4.6|) : 



In analogy to (5.9b) we would then conjecture that the r-th order microphysical decaying 
state is described by the state operator 



n=0 ^ ^ ■ k=0 ^ ^ 



r-1 , 

27rryf ^ 

n=0 ^ 



(up to a normalization factor which will have to be determined by normalizing the overall 
probability to 1). Since (5.9b) is postulated to be the zeroth order Gamow state representing 
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a resonance, ( |5.10| b) is conjectured to be the r-th order Gamow state.0 

Whether the microphysical state of the (hypothetical) quasistationary microphysical sys- 
tem is always represented by the mathematical object ( |5.10| b) or whether also each individual 

k=0 ^ ^ k=0 

n = 0,l,...,r-l , (5.11b) 

has a separate physical meaning, is not clear. So far it is not even certain that higher order 
poles of the S-matrix describe anything in nature (though there are no theoretical reasons 
that exclude these isolated singularities of the S-matrix.) But if these hypothetical objects 
do exist, the r-th order pole is associated with a mixed state (|5.1CI| b) whose irreducible 
components are given by ( [5.11| b). E.g., for the case r = 2 (second order pole at z^) we have: 



and 



and 



iy(o) = |^^)(o){o)(-z^| (5.12) 



Ty« = r(|z^)W + W(-z^l) (5.13) 



W = 27rr 



{\z-r ^'H-zn\ - 2.r(|.-)(°) «(-.^J + |.-)« W(-.^l)) (5.14) 



This means that the conjectural physical state associated with the r-th order pole is a mixed 
state W, all of whose components W^"'\ except for the zeroth component W^^\ cannot be 
reduced further into "pure" states given by dyadic products like \z'^)^''^ zr\. This is 
quite consistent with our earlier remark that the label k is not a quantum number connected 
with an observable (like the suppressed labels 62, • • • ,bn)- Therefore a "pure state" with a 
definite value of k, like l-z^)*-'^^ ^'''>{~zr\, k > 1, does not make sense physically. A physical 
interpretation could only be given to the whole r-dimensional space AizR, (|4.14]) . The 
individual W^^\ n = 0, 1, 2, ...,r — 1, act in the subspaces A^i*^^ C which are 

spanned by Gamow vectors of order 0, 1, . . . , n (Jordan vectors of degree n + 1, i.e. {H^ — 
zr)^~^^A4^J^j^ = 0). There the question is, whether there could be a physical meaning to 
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We want to mention that mathematically there is an important difference between ( ^.10^ ) and ( ^.f O b) 



because (■0 \z )('^) (" '^)(+2;|0+) are analytic functions for z in the lower half-plane, whereas the 
ii'^lz-)^''') are not. 
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each W^""^ separately, or whether only the particular mixture W given by ( 5.1(J| b) can occur 
physically. 

Though the quantities l^;^^)*-'^'' ^''H~zr\ will have no physical meaning, even if higher or- 
der poles exist, they have been considered |jl9[ and their time evolution is calculated in a 
straightforward way from (|4.35|) : 

e-''''''\z-.y^'^^'^^-zn\e''''= (5.15) 

1 1 
/! ml 



k k 



1=0 m=0 

This time dependence (as well as the time dependence in ( ^4.35| )) is reminiscent of eq. ( |4.9| ) 
in the reference of M. L. Goldberger and K. M. Watson 0. 

It shows the additional polynomial time dependence, that has always been considered 
an obstacle to the use of higher order poles for quasistationary states. A polynomial time 
dependence of this magnitude (of the order of r = ^) should have shown up in many 
experiments. 

We now derive the time evolution of the microphysical state operator defined in (|5.11[ b) 
using the time evolution obtained for the Gamow- Jordan vector in ( [4.36[ ). It will turn out 
that this operator, whose form was conjectured in analogy to the pole term ( |5.11| a), will have 
a purely exponential time evolution. This was quite unexpected. 

Inserting (|4.36| a) and ( [4.36| b) into 

l^W(t) = ^-^H^t^in)^^m = I! ^ Q c"^^^ * | ) Zr\ e^^* (5.16) 
we calculate: 

fc=0 /=0 m=0 V / V / V / 
n n—m n—m / \ / j \ / j\ 

EET.il) ii) [m ) i-^tr'i^tr-'-"'\^Rf' zr\ (5.17) 

m—n ;— n ; — ; \ / \ / \ / 



n! 

m=0 /=0 k=l 

n n—m n—m 

-rt 



■ m=0 Z=0 fc=« ^ ^ ^ ^ ^ ^ 

^ E (:) E h '") '-'(-^hI E i'T-i') i-nr-"- 

m=0 ^ ^ i=0 ^ ^ k=l ^ ' 
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k 




Figure 3: For the case n = 3, the summation terms labeled by the parameters k, /, and m 
are displayed as dots in the diagram to show that the summations of lines 2 and 3 of ( p. IT] ) 
both contain the same terms. 

In going from the second to the third line, the order of summation has been changed, by 
keeping the same terms, as displayed in fig. 3 for the case n = 3. In going from the third to 
the fourth line one uses the identity 

fn\ fk\ fn — k\ f n\ fn — m\ fn — m — l\ io\ 
[k)\l)\m ) = [m) \ I )[ k-l ) (^-^^^ 

where ( ? ) = kUn-kV. binomial coefficients. Since the indices labeling the Gamow- Jordan 
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vectors do not depend upon k, the sum over k may be performed using the binomial formula: 



Efn — m — t\, .^._i, .n-k-m (■. ■,\n-m-l Jl ioT I — H — TTl \ ^ 
[ k-l ) ^-'^^ =(.t-^t) =|o for/^n-m/ = '^'" 

(5.19) 

Inserting ( [5.19[ ) into the fourth line of ( |5.17| ) and performing the sum over / then gives: 



■pn " / \ 

= e"^*— 5Z (^) = e-^*W^^"nO) ; t>0 (5.20) 

This means that the complicated non-reducible (i.e. "mixed") microphysical state oper- 
ator defined by ( |5.11| b) has a simple purely exponential semigroup time evolution, like 
the zeroth order Gamow state (5.9b). This operator is probably the only operator formed by 
the dyadic products l^^)^™-* ^'-'("-z/?! with m, 1 = 0, 1, ■ ■ ■ , n, which has a purely exponential 



time evolution. Thus W^^'' of eq. (|5.11[ b) is distinguished from all other operators in Ai^J^^ 



The microphysical decaying state operator associated with the r-th order pole of the 
unitary S-matrix is according to its definition ( |5.10| b) a sum of the W^'^\ Because of the 
simple form (|5.2CI|) (independence of the time evolution of n) this sum has again a simple 
and exponential time evolution 

Wit) = e-'^"'We'^' = e-^'W ; t > . (5.21) 

Thus we have seen that the state operator which we conjecture from the r-th order pole 
term describes a non-reducible "mixed" microphysical decaying state which obeys an exact 
exponential decay law. 

We can return to the question that we started with when we set out to conjecture the 
state operator for the (hypothetical) microphysical state associated with the r-th order S- 
matrix pole: What is the probability to register at the time t the decay products \ (or 

in general A = |V'j~)(V'j~l) if at t = the microphysical state was given by W of ( ^.lOj b)? 



From (|5.21| ) we obtain 

PaH) = TT{AW{t)) = e-^*Tr(AH^) = e-^*PA(0) (5.22) 

or in the special case of A = 

P^-{t) = (niy(t)l^-) = e-^'{^~\W\^-) (5.23) 

This is exactly the same result as the result ( p.3| b) for the microphysical state W'-^ asso- 
ciated with the first order pole of the S-matrix and the result which is in agreement with 
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the experiments on the decay of quasistationary states. It is, however, important to note 
that in our derivation of ( ^.20|) and (|5.23|) we proceeded in a very specific order. We first 
derived ( |5.2CI| ) from ( |4.36| a) and ( |4.36| b) and then calculated the matrix elements with ip' 
and not vice versa in order to avoid problems with the analyticity. 



6 Conclusion 

Vectors that possess all the properties that one needs in order to describe a pure state of 
a resonance have been known for two decades. These Gamow vectors tp'^ are eigenvectors 
of a self-adjoint Hamiltonian with complex eigenvalues Eji — iV /2 (energy and width), they 
are associated with resonance poles of the S-matrix, they evolve exponentially in time, and 
they have a Breit-Wigner energy distribution. They also obey an exact Golden Rule, which 
becomes the standard Golden Rule in the limit of the Born approximation. The existence 
of these vectors in the rigged Hilbert space allows us to interpret exponentially decaying 
resonances as autonomous microphysical systems, which one cannot do in standard Hilbert 
space quantum mechanics. 

The mathematical procedure by which these Gamow vectors had been introduced suggests 
a straightforward generalization to higher order Gamow vectors which are derived from higher 
order S-matrix poles. We have shown in this paper that the r-th order pole of a unitary 
S-matrix leads to r generalized eigenvectors of order /c = 0, I,-- - , r — 1. These fc-th order 
Gamow vectors are Jordan vectors of degree {k -\- 1) with complex eigenvalue — iV /2. 
They are basis elements of a generalized eigenvector expansion. But their time evolution has 
in addition to the exponential time dependence also a polynomial time dependence, which 
is excluded experimentally. However, the generalized eigenvector expansion suggests the 
definition of a state operator for microphysical decaying states of higher order. These state 
operators cannot be expressed as dyadic products of generalized vectors. But these state 
operators have a purely exponential time evolution. 



There has been a lot of interest in the Jordan blocks for various applications (see e.g. ||TT 



T2| , |T3|, |I5|, |T9], |16]). Here it has been shown that Jordan blocks arise naturally from higher 
order S-matrix poles and represent a self-adjoint Hamiltonian by a complex matrix in a 
finite dimensional subspace contained in the rigged Hilbert space. Although higher order 
S-matrix poles are not excluded theoretically, there has been so far very little experimental 
evidence for their existence, because they were always believed to have polynomial time 
dependence. Since we have shown here that their non-reducible state operator evolves purely 
exponential in time there is reason to hope that these mathematically beautiful objects will 
have some application in physics. 
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